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ABSTRACT 

We suggest that the abrupt switch, from hierarchical clustering on scales ~ 0.04 pc, to bi- 
nary (and occasionally higher multiple) systems on smaller scales, which Larson has deduced 
from his analysis of the grouping of pre-Main-Sequence stars in Taurus, arises because pre- 
protostellar gas becomes thermally coupled to dust at sufficiently high densities. The resulting 
change from gas-cooling by molecular lines at low densities to gas-cooling by dust at high 
densities enables the matter to radiate much more efficiently, and hence to undergo dynamical 
fragmentation. 

We derive the domain where gas-cooling by dust facilitates dynamical fragmentation. 
Low-mass (~ Mq) clumps - those supported mainly by thermal pressure - can probably ac- 
cess this domain spontaneously, albeit rather quasistatically, provided they exist in a region 
where external perturbations are few and far between. More massive clumps probably require 
an impulsive external perturbation, for instance a supersonic collision with another clump, 
in order for the gas to reach sufficiently high density to couple thermally to the dust. Im- 
pulsive external perturbations should promote fragmentation, by generating highly non-linear 
substructures which can then be amplified by gravity during the subsequent collapse. 

Key words: ISM: clouds - dust processes - Stars: formation 



1 INTRODUCTION 

Dynamical fragmentation (i.e. the growth of substructure during 
near-freefall collapse) is widely believed to play an important role 
in star formation, both as a means of forming star clusters, and 
as a means of forming binary systems (e.g. Bate, Bonnell & Price 
1995; Bodenheimer, 1978; Bodenheimer & Boss, 1981; Bonnell 
& Bastian 1993; Bonnell & Bate 1994; Boss, 1991, 1993; Boss 
& Bodenheimer 1979; Burkert & Bodenheimer 1996; Chapman 
et al., 1992; Hoyle, 1953; Larson, 1972; Miyama, 1992; Miyama, 
Hayashi & Narita, 1984; Monaghan, 1994; Monaghan & Lattanzio, 
1991; Myhill & Kaula, 1992; Nelson & Papaloizou, 1993; Pringle, 
1989; Sigalotti & Klapp, 1994; Turner et al., 1995; Whitworth et 
al., 1995). 

Recently Larson (1995) has analyzed the grouping of pre- 
Main-Sequence (pre-MS) stars in Taurus, and identified a critical 
length-scale L cril ~ 0.04 pc. On larger scales the pre-MS stars 
appear to be clustered in a self-similar hierarchy, with the mean 
surface-density of neighbours N decreasing with angular separa- 
tion 8 as TV oc 9~ 62 . On smaller scales the pre-MS stars have 
on average just one companion, which is usually tightly bound in 
a binary system; here the mean surface density of companions de- 
creases as N oc (9 -215 . Larson suggests that these binary (and 
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occasionally higher multiple) systems are formed by dynamical 
fragmentation. He proposes that L aii is simply the Jeans length 
appropriate for the conditions prevailing in Taurus, and that the 
corresponding Jeans mass is ~ Mq. He then points out that if 
there is a preponderance of marginally Jeans-unstable clumps with 
M ~ Mq collapsing and fragmenting into binary and triple sys- 
tems, this will yield protostellar masses in the range 0.3 to 0.5 Mq, 
i.e. very close to the characteristic mass of the Initial Mass Func- 
tion, ~ 0.4M©. 

However, there is a large range of densities in Taurus, and 
presumably an even large range of conditions in the various star- 
formation regions which form the entire ensemble of stars defining 
the Initial Mass Function. Therefore there must be a wide range of 
Jeans lengths. 

Here we argue that the onset of fragmentation at the scale L ail 
is in fact due to the density becoming just high enough for the gas 
to couple thermally with the dust, and thereby greatly increase its 
cooling efficiency, by availing itself of the dust's continuum opac- 
ity. This circumstance defines a length-scale ~ L a - a and a mass- 
scale ~ Mq which are only weakly dependent on metallicity and 
environmental factors. 

We concur with Larson's inference that dynamical fragmen- 
tation will then tend to deliver protostars with masses around the 
peak in the IMF at O.4M . 

Dynamical fragmentation can only occur if the gas cools ef- 
ficiently by radiating. Specifically, the gas must be able to radiate 
away - on a dynamical timescale or faster - the internal energy gen- 
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erated by gravitational compression (Hoyle 1953, Rees 1976, Low 
& Lynden-Bell 1976, Silk 1977). At high densities this requires 
continuum opacity, and for cold (T ~ 10 2 K, say) gas having fi- 
nite metallicity this opacity is supplied by dust. 

Recently, Larson's analysis of the grouping of pre-MS stars 
in Taurus has been confirmed by Simon (1997). Simon has also 
repeated this analysis for the p Ophiuchus and Orion-Trapezium 
star-formation regions. He finds that the surface-densities of neigh- 
bours and companions obey similar power laws, TV cc g-°- 5±0 - 2 



on large scales, and N oc 9 



on small scales. However, 



he finds that the switch occurs at different projected separations: ~ 
0.06 pc in Taurus, ~ 0.025 pc in p Ophiuchus, and ~ 0.002 pc in 
Orion-Trapezium . 

Simon argues that this variation is difficult to attribute to sys- 
tematic changes in the Jeans length between the different star- 
formation regions. We propose that some of this variation arises be- 
cause the different star-formation regions have very different large- 
scale surface densities. Proceeding from Taurus, through p Ophi- 
uchus, to Orion-Trapezium, the large-scale surface-density of stars 
increases significantly, and so the switch must shift to smaller pro- 
jected separations, even if the intrinsic distribution of binary sepa- 
rations is identical. 

Kitsionas, Gladwin and Whitworth (1998) have also re- 
analyzed Taurus, with a larger sample of pre-MS stars. They have 
repeated this analysis for the Chamaeleon I star-formation region. 
Power-law fits to N(6) give exponents similar to those obtained by 
Larson (1995) and Simon (1997). For Chamaeleon I the switch is 
at ~ 0.03 pc, as would be expected on the basis of its high overall 
surface-density. 

In Section 2 and Appendix A we derive an expression for the 
compressional heating rate, TL comp , in a recently formed fragment. 
In Section 3 we evaluate thermal coupling between gas and dust. In 
Section 4 and Appendix B we derive an expression for the cooling 
rate, £ d ust, due to thermal emission from dust. Appendix C demon- 
strates that for fragments having radius R ^ 0.1 pc, gas-cooling by 
dust dominates over gas-cooling by lines. In Section 5 we derive 
an expression for the mean (frequency-averaged) optical depth due 
to dust. In Section 6 we combine the results from the previous sec- 
tions to define the domains on the (R, M) and (p, T) planes where 
gas-cooling by dust facilitates dynamical fragmentation. In Section 
7 we discuss the results, and in Section 8 we summarize our main 
conclusions. 
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In Appendix A we develop a simple time-dependent model of Jeans 
instability which enables us to estimate J (~ 2.36), and M(S). 
For the purpose of numerical estimates, we set S = 0.2 (up to 5 
fragments). 

For completeness, we also include a cosmic-ray heating rate 



Hcr = r M 



(3) 



with r„ ~ 2 x 10~ 4 erg s _1 g^ 1 . This corresponds to a cosmic- 
ray ionization rate X C r — 6 x 10 -17 s _1 . Eqn. ^ makes no al- 
lowance for the attenuation of cosmic rays in the interior of a dense 
fragment. However, this can only reduce the cosmic-ray heating 
rate, and it turns out that even the unattenuated rate is not very im- 
portant in the regime with which we are concerned, i.e. the regime 
where gas-cooling by dust facilitates dynamical fragmentation. 



3 THERMAL COUPLING BETWEEN GAS AND DUST 

If the gas in a fragment is thermally coupled to the dust, the contin- 
uum opacity of the dust affords the fragment an efficient means of 
radiating away the internal energy generated by gravitational com- 
pression. Here we formulate the rate C at which the gas transfers 
thermal energy to the dust. 

If a fraction Z im by mass of the interstellar medium is in the 
form of spherical dust grains having radius r dusl and internal density 
p dllsl , the rate of transfer of thermal energy to the dust is given by 



Cm = 



1/2 



3pZ dl 



M /8fcT\ J 
rh \ urn ) 4r d „ st pdu S t 



exp 



-75K 



T 



k(T - r dusl 

(7^1) 



(4) 



For the purpose of numerical estimates we put r dlls i 



10" 



and p dusl = 3 g cm . The term in braces, {}, represents the ac- 
commodation coefficient, i.e. the mean inelasticity of collisions be- 
tween gas particles and dust grains (Hollenbach & McKee 1979). 
For the low temperatures (T ~ 100 K) with which we are mainly 
concerned here, 7 ~ 5/3, since the rotational degrees of freedom 
of H2 are not excited. 



2 COMPRESSIONAL HEATING 

For an approximately spherical fragment which has just become 
Jeans unstable, the mass M,, radius R,, temperature T and mean 
gas-particle mass rh are related by 

GM, JkT 2 

—5— = — = Ja„ , (1) 
R s m 

where J is a numerical factor of order unity (see below), k is Boltz- 
mann's constant, and a a = [kT/rh] 1 ' 2 is the isothermal sound 
speed. 

If the gas remains isothermal, then by the time this frag- 
ment has contracted sufficiently for sub-fragments of mass SM S 
to become Jeans unstable, the original fragment has radius R = 
S 2 ^ 3 Ri, and is collapsing at speed dR/dt = — M(S)a . Conse- 
quently the compressional heating rate is 

dt dt 



4 RADIATIVE COOLING 

The maximum radiative cooling rate is obtained by assuming that 
the fragment radiates like a blackbody, i.e. 

£bb = 4 7ri? 2 a SB T 4 . (5) 

This is a maximum cooling rate on four counts. 

First, it assumes that the collapsing fragment exists in isola- 
tion. In reality, part of the luminosity from the fragment will be 
cancelled by the flux of background radiation incident on the sur- 
face of the fragment. The background radiation field will consist of 
a cosmic component at T ~ 3 K, plus a contribution from local 
sources such as stars, protostars, circumstellar cocoons and discs, 
warm molecular clouds. A detailed evaluation of this contribution 
would be extremely model-dependent. For the purpose of this pa- 
per, we represent the background with a single blackbody compo- 
nent at temperature T„ ~ 3 K. 
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Second, it assumes that the cloud has continuum opacity at all 
wavelengths - or, at the very least, across a broad range of wave- 
lengths near the peak of the blackbody spectrum. 

Third, it assumes that, at these wavelengths around the peak, 
scattering is negligible compared with absorption and emission. 

Fourth, it assumes that, at these wavelengths around the peak, 
the cloud is neither optically thin, nor very optically thick. 

If we define a mean optical depth r(T) between the centre and 
the surface of the fragment (i.e. an optical depth averaged over the 
frequencies which contribute to radiative cooling of the fragment 
at temperature T), the radiative cooling rate is given approximately 
by 



4 TV R (7 SB 



T 4 f (T)r c „ 



(r(T) + l)(f(T) + r crll ) 



T 4 f (T )r cril 



(f(T ) + l)(f(T u )+r tr „) 



(6) 



Here, r tri , is a critical optical depth, typically ~ 20, and the second 
term in braces simply ensures that the cloud must be hotter than the 
background to cool radiatively. Eqn. (^) is constructed so as to give 
the correct asymptotic forms. For example, if T„ = 0, then in the 
optically thin limit, 



— ► 4nR A <T S &T' 
at modest optical depths, 

C rili > 47ri? 2 a SB T 4 

and at high optical depths, 

C aL , — > 47r7? 2 crsB 



-(T) (t(T) < 1) 



(1 < f (T) < r„,) 



f(T) 



(r(T) > r„„) 



This last form is derived in Appendix B, and takes account of the 
fact that high luminosities can only be transported through the in- 
terior of an optically thick fragment if there is a sufficiently large 
temperature difference between the centre of the fragment and its 
surface. To take account of this, we require that no more than a 
fraction / of the fragment's mass be hotter than FT, where T is 
now specifically the surface temperature. Hence the precise value 
of Ten depends on our choice of / and F, and also on f3 (the 
emissivity index): r tri , = r cr ,,(/ ', F, (3). Substituting / = 1/3, 
F = 2 (i.e. at least 2/3 of the fragment is below 2T), we obtain 
Tcrii = 47.9 & 30.8 for j3 = 1 & 2 respectively (see Appendix 
B). For the purpose of numerical estimates we shall adopt /3 = 1.5 
and Tent = 38.3. 



5 DUST OPACITY 

The mean dust optical depth between the centre and the surface of 
a spherical, uniform-density fragment is 



F (T dllsl ) = 



9MZ^Q(T^) 
167r_R 2 r du . 



An* 



where Q(Tdiisi) is the Planck-mean absorption efficiency of a dust 
grain. Since our main concern here is with long wavelength radia- 
tion from dust at low temperatures, we adopt the prescription 



Q(T dust ) 



0.0004 



r llu 



10K 



1 



This corresponds to a mean mass-opacity coefficient K,(T iaa i 
3Z du!t g(T d „ s ,)/4r du!t p d „ 8t ~ 10 cm 2 g- 1 Z duA [T dllsl /10K]' 3 , 
hence a mean optical depth 



and 



2.4 cm' 



z M 



10K 



(7) 



(cf. Hildebrand (1983), who recommends Q(X — 250/im) = 
0.0004 and n(\ = 250pm) = 0.1 cm 2 g^ 1 ). 

To account for dust evaporation, we multiply Eqn. (Q) by a 
factor 



0.4 exp 



100K 



T d „ 



0.6 exp 



1000K 



T,\,, 



The second term represents evaporation of the volatile grain mantle 
at temperatures of order 100 K, and the third term represents evap- 
oration of the refractory core at temperatures of order 1000 K. This 
rather ad hoc prescription is adequate for our purposes, since we 
are mainly concerned with gas-cooling by dust in the regime where 
the dust is cool and evaporation is negligible. 



6 GAS-COOLING BY DUST DURING DYNAMICAL 
FRAGMENTATION 

The domain in which gas-cooling by dust facilitates quasi-freefall 
collapse and dynamical fragmentation (hereafter the dust-domain) 
is given byPl 



He 



+ Hctt ~ {Cg-d + -Cad 1 } 



(8) 



The dust-domain can be defined on both the (R, M) and the (p, T) 
planes, using Eqn. (|l|) and p = 3M/AtvR s to convert from one to 
the other. 

In order to simplify the analysis, we put T dus , = T = 10 K 
in C„_ d (Eqn. (Q)), and T dusl = T in £„, d (Eqn. To justify this 
we note that the section of the dust-domain boundary where C g d is 
important is quite distinct from the section where £ Dld is important, 

viz. 

Consider first the large-ii limit of the dust-domain, where 
(L_ d Cai , so Cg-d dominates the righthand side of Inequality 
<jj8J>, i.e. the critical process is the transfer of thermal energy from 
the gas to the dust. Under this circumstance the thermal balance 
of the dust involves a much more rapid energy turnover than that 
of the gas (e.g. Whitworth & Clarke 1997). Hence the heat input 
from thermal coupling to the gas makes a negligible contribution 
to the thermal balance of the dust. The dust temperature T dusl is de- 
termined by the background radiation field, and adopts a value of 
order the background temperature T . 

Next consider the small-i? limit of the dust-domain, where 
C s - d 2> Cms , so Cms dominates the righthand side of Inequality 
<j|s|> , i.e. the critical process is the emission of radiation by the dust 
and its escape from the fragment. At these high densities the gas 
and dust are closely coupled. A very small temperature difference 
suffices to transfer the compressional heating of the gas to the dust, 
and we can approximate T dus , ~ T. 

In between these two limits, the dust-domain boundary is de- 
termined at high temperatures by dust evaporation, and at low tem- 
peratures by the fact that the gas cannot cool below the background 
temperature. Consequently there is no need to formulate accurately 
the circumstances under which C g d ~ £ rad . 

Figure 1 shows contours of constant 



{Hcomf + Hcr) (CL/ + C r .J) 



° The two heating processes operate simultaneously, and so their rates add 
in parallel. The two cooling processes operate sequentially, and so their 
rates add in series. 
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Figure 1. The (R, M) plane. The outer contour is the limit of the dust- 
domain, where gas-cooling by dust enables the gas to fragment dynamically. 
The other contours show where the ratio of compressional heating to gas- 
cooling by dust equals 0.32, 0.10 & 0.03. The letters refer to the different 
processes limiting the dust-domain, as discussed in the text. The line of open 
circles represents the average observed properties of clumps. The large dot 
marks the point (0.04 pc,Mq), i.e. the critical entity identified by Larson. 



Figure 2. The (p, T) plane. The outer contour is the limit of the dust- 
domain, where gas-cooling by dust enables the gas to fragment dynami- 
cally. The other contours show where the ratio of compressional heating to 
gas-cooling by dust equals 0.32, 0. 10 & 0.03. The line of open circles repre- 
sents the average observed properties of clumps. The dashed line represents 
the background temperature. 



on the (R, M) plane, for f3 = 1.5, T„ = 10K, and Z, us , = 0.01. 
The region inside the outer contour is the dust-domain. Contours 
are drawn for 1Z = 1.00, 0.32, 0.10 & 0.03, to emphasize how 
efficient gas-cooling by dust becomes inside the dust-domain. 

Starting at the righthand side, and working round in a clock- 
wise sense, the dust-domain limits are caused by (A) gas-dust cou- 
pling in the high-temperature (elastic) limit; (B) gas-dust coupling 
in the low-temperature (perfectly inelastic) limit; (C) the back- 
ground temperature acting as a minimum temperature (T > T„); 
(D) the dust being optically thin; (E) the dust optical depth being 
intermediate (1 ~ f ~ r cril ); (F) the dust being optically thick; 
(G) dust evaporation. Partial dust evaporation actually promotes 
fragmentation in the high density limit by reducing the dust op- 
tical depth, but it inhibits gas-dust coupling; this is why the dust- 
domain boundaries curve towards higher densities as they approach 
the evaporation limit. The same contours are presented on the (p, T) 
plane in Figure 2. The dashed line indicates the background tem- 
perature T„. 

The most relevant part of the dust-domain is probably the part 
corresponding to low temperatures, T ^50K, but we have plotted 
the whole domain in order to avoid the extremely model-dependent 
complexities of predicting the gas temperature precisely. 

Of course, to be important in this context gas-cooling by dust 
must not only be able to handle the net heating of a fragment, it 
must also be able to handle it more effectively than, or at least as 
effectively as, gas-cooling by lines. Figure 1 shows that the dust- 
domain has R JSO. lpc, and Appendix C demonstrates that for pro- 
tostellar fragments with R SS>0. lpc, gas-cooling by dust is more 
effective than gas-cooling by lines. 



7 DISCUSSION 

On Figure 1, we have also plotted the mean observed distribution 
of clump masses and radii, as collated by Chieze (1987). This dis- 
tribution approximates to 



A I 



M G 



R 



O.lpc 



and is marked by a line of open circles on Figure 1 . If the clumps 
were thermally virialized, we could use Eqn. ([!]) to convert this 
distribution into 



T 



8.5 K 



10" 20 g cm" 3 



and this is marked by the line of open circles on Figures 2. 

However, all but the very low-mass clumps are supported by 
non-thermal motions, with the gas being at rather low temperatures, 
typically in the range 10 to 15 K, and seldom above 50K. Thus, 
although the locus of clumps in Figure 1 appears to run rather close 
to the boundary of the dust-domain, it is only the low-mass clumps 
- those supported mainly by thermal pressure - that are really close 
to the dust-domain. Even these may be unable to access the dust- 
domain without some external perturbation. 

The gas in a high-mass clump can probably only enter the 
dust-domain if it is compressed and heated by a shock, for exam- 
ple following a collision with another clump of comparable mass, 
or as a consequence of being overrun by an expanding nebula (HII 
region, stellar-wind bubble or supernova remnant). Gas-cooling by 
lines will play an important role in such shocks (as a means of 
removing the dissipated energy), but Appendix C shows that the 
shock-compressed gas will then only be able to condense further 
and fragment if it is dense enough to couple thermally to the dust. 
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A number of observations might be explained on the basis of 
the dust-domain and its relation to the locus of observed clumps. 



7.1 Low-mass clumps 

Low-mass (~ Mq) clumps are supported mainly by thermal pres- 
sure (Myers 1983). Therefore, in the absence of external perturba- 
tions, they should relax to hydrostatic equilibrium. However, Fig- 
ures 1 & 2 suggest that such clumps need a relatively small external 
perturbation - and possibly no perturbation at all - to enter the dust- 
domain. Consequently we might expect a large fraction of stars to 
be formed from the collapse of such clumps. If dynamical fragmen- 
tation leads typically to the formation of a binary or triple system, 
with most of the clump mass ending up in the component stars, 
then the Initial Mass Function (IMF) should peak in the range 0.2 
to 0.5 Mq, as it apparently does (Scalo 1986). This last point was 
emphasized by Larson (1995). 

Moreover, if an external perturbation is involved, this should 
greatly increase the chances of permanent fragmentation during the 
ensuing collapse. This is because the clump would be delivered into 
the dust-domain with highly non-linear substructure, and gravity 
would then rapidly amplify this non-linear substructure during the 
subsequent quasi-freefall collapse. 

By contrast, if it is possible for a clump to condense quasistat- 
ically into the dust-domain, it arrives there in a rather too well orga- 
nized state, and collapse seems likely to lead to a single star. Under 
this circumstance, binary formation may require some sort of in- 
trinsic instability, for example in a circumstellar disc. 

Support for the notion that even low-mass clumps either re- 
quire an external perturbation to trigger their collapse, or con- 
dense rather quasistatically into the dust-domain before collapsing 
to form stars, comes from the observation (Myers, private commu- 
nication) that the majority of starless low-mass clumps are located 
in the outer reaches of Taurus, where presumably external pertur- 
bations are few and far between. 

Further support comes from the observation that in low-mass 
clumps associated with IRAS sources, the IRAS source (i.e. the pu- 
tative newly-formed star) is normally significantly displaced from 
the centre of the (presumably placental) clump (Beichman et al. 
1986, Clark 1987). 



7.2 High-mass clumps 

High-mass (~ IOMq) clumps are supported mainly by non- 
thermal motions. From Figure 2, it appears that they are unable 
to enter the dust-domain, without a substantial external impulse - 
most likely a shock, due to a collision with another comparable 
clump, or due to being overrun by an expanding nebula. A suffi- 
ciently strong shock would compress the gas to the point where it 
became thermally coupled to the dust. If our hypothesis is correct, 
this would lead to efficient fragmentation and the formation of a 
star-cluster. 

There is evidence that molecular clouds have an hierarchical 
structure, i.e. clumps nested within clumps. If this is so, a high- 
mass clump should be envisaged as an ensemble of lower-mass 
clumps. Consequently, even when star formation is triggered in 
a high-mass clump (by an external perturbation), most of the gas 
ending up in stars will probably reside in low-mass (~ Mq) sub- 
clumps, and this should lead to a majority of the stars formed hav- 
ing masses in the range 0.2 to 0.5 Mq. 




-4 -2 2 

Log(R/pc) 

Figure 3. The (R, M) plane. The effect on the dust-domain of varying the 
far-infrared emissivity index /3. The heavy contour delimits the dust-domain 
when (3 = 1. The light contour delimits the dust-domain when /3 = 2. Note 
that the two contours only diverge at small R. 

7.3 M miB and M mx 

Several authors (Rees 1976, Low & Lynden-Bell 1976, Silk 1977) 
have determined the minimum mass for star formation M m[B on the 
basis of opacity-limited fragmentation. With reference to Figure 1, 
this is equivalent to identifying M min with the smallest mass falling 
inside the dust-domain, and it yields a value M mill ~ O.OIMq. 
However, to reach this part of the dust-domain from the conditions 
obtaining inside a clump would require either (i) several consecu- 
tive fragmentation episodes, or (ii) a single fragmentation episode 
generating a very large number of pieces. Numerical simulations 
(Larson 1978) suggest that dynamical fragmentation is a one-stage 
process, i.e. that during freefall gravity only amplifies the substruc- 
ture which is present at the outset. This implies that (i) cannot 
work, and that (ii) requires an external perturbation with improb- 
ably high power on high spatial frequencies. Therefore we do not 
expect many stars with masses ~ O.OIM©. 

The largest mass falling inside the dust-domain is M m!lx ~ 
1OOM . We speculate that this is because clumps having higher 
mass cannot (coherently) enter the dust-domain, no matter what ex- 
ternal perturbation they experience, and therefore cannot collapse 
to form stars. 



7.4 Dependence on /3, T;, and Z dusl 

In Figure 3 we show how the dust-domain is altered on the (R, M) 
plane, if one adopts a different far-infrared emissivity index for the 
dust, i.e. /3 = 1 or f3 = 2. The only significant change is on the 
small-i? boundary of the dust-domain, where the dust is optically 
thick to its own radiation. Since we are here more concerned with 
the large-i? boundary, this variation will not affect our conclusions 
significantly. 

In Figure 4 we show how the dust-domain is altered on the 
(R, M) plane, if one adopts a different background temperature, 



© 0000 RAS, MNRAS 000, 000-000 



6 A. P. Whitworth, H. M. J. Boffin and N. Francis 



2 - 



o 




-4 



-2 



Log(R/pc) 



Figure 4. The (R, M) plane. The effect on the dust-domain of reducing the 
metallicity to Z dusl = 0.003, heavy contour; and reducing the background 
temperature to T = 5 K, light contour. 



i.e. T = 5 K (instead of T = 10 K), or a different metallicity, i.e. 
Zfa St — 0.003 (instead of Z dusl = 0.01. Changing the background 
temperature has little effect on the boundary of the dust-domain 
where it is close to the locus of clumps. 

Reducing the metallicity moves the dust-domain to smaller R- 
values. Unless the locus of clumps is also different, in the sense of 
being characterized by a larger mean column-density, this increases 
the gap between the dust-domain and the locus of clumps. Conse- 
quently, where the metallicity is well below solar, it should be much 
harder for star formation to occur spontaneously. This might mean 
that star formation is more dependent on an external impulse, and 
hence that star formation is more likely to occur in bursts. This 
could dictate the pattern of star formation in a galaxy having low 
metallicity, for example a young galaxy or an irregular galaxy - in 
the sense of making star formation more burst-like. 



8 CONCLUSIONS 

We have evaluated the conditions under which gas-cooling by dust 
enables a clump to collapse and fragment dynamically. We describe 
these conditions as the dust-domain. We have also shown that in 
the dust-domain gas-cooling by dust dominates over gas-cooling 
by dust. 

The observed properties of clumps in star-forming molecular 
clouds are such that low-mass (~ M©) clumps are close to the 
dust-domain. Consequently only a small external perturbation is 
needed to propel them into the dust-domain - and possibly no per- 
turbation at all. Therefore, if such clumps fragment into binary or 
triple systems, with most of their mass ending up in the components 
of these systems, they will spawn stars having masses concentrated 
in the range 0.2 to 0.5 M®, which is the peak of the IMF (cf. Larson 
1995). 

In isolated quiescent regions, low-mass cores may survive for 
long periods, for lack of a sufficient external impulse to trigger 



(or accelerate) their collapse. They would be observed as starless 
clumps. 

High-mass clumps require a substantial external impulse to 
reach the dust-domain and start collapsing. Collisions with other 
clumps and/or shocks due to expanding nebulae can provide a suffi- 
cient impulse. This may be why star formation is (apparently) more 
efficient in dynamically agitated regions. 

The dust-domain may also define the minimum and maximum 
masses for star formation, M min ~ 0.01M Q and M max ~ 1OOM , 
but the majority of stars should fall well within these limits. 

In regions of low metallicity, for instance irregular galaxies, 
star formation may be more prone to occur in bursts. 



Appendix A. Isothermal fragmentation 

Consider a uniform medium with density p and isothermal sound 
speed a . Now suppose that a spherical fragment having radius R , 
and hence mass M = 47ri? s p /3, attempts to condense out of 
the background by contracting homologously, so that it remains 
spherical with instantaneous radius R, and uniform with instanta- 
neous density p — p„ (R/R ) 3 . Further suppose that the isother- 
mal sound speed is constant, a = a„, and that the external medium 
exerts a constant pressure P cx , = p a 2 at the boundary of the frag- 
ment. Radial excursions of the fragment are controlled by a poten- 
tial function 4>(R), in the sense that d 2 R/dt 2 = - d(j>/dR and 
(dR/dt) = [~2(j)(R)] 1/2 . 4>(R) is derived (Whitworth 1981) 
from the energy equation 



(9) 



where K, = (3/W)M„(dR/dt) 2 is the radial kinetic energy, Q — 
— 3GM 2 /5R is the self-gravitational potential energy, and P iM = 
3M a 2 /4irR :i is the internal pressure. 4>(R) can be written in the 
form 



R) _ _15_ (FL^ 2 (J^y 1 _ 



R_ 

,Ro 



(10) 



where the first term on the righthand side represents self-gravity, 
the second term internal pressure, the third term external pressure, 
and 



a = 



3 

2 7/3 



TvGp„ 



1/2 



(ID 



Stable equilibrium states (d(j>/dR = 0, d 2 <f>/dR 2 > 0) exist only 
for 7? < R s ; for R > R, the fragment contracts indefinitely. 
Therefore R s can be identified as the leans radius. The leans mass 
is given by 



AnRfp 



3 2 tt / M 3/2 
2 5 UgJ p i/2 ' 



and so 



J = 



GM, 
R,a 2 a 



15 

2 8 / 3 



2.3623 



(12) 



(13) 



If we consider a marginally leans-unstable fragment with 
mass Mj, starting to condense out from radius R s , then by the time 
sub-fragments with mass SM S become leans unstable, the radius 
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of the initial fragment has decreased to S 2 ^ 3 R S . From Eqn. (|Io[), it 
follows that the Mach Number of the contraction is given by 



M(S) 



dR/dt 



(14) 



Representative values of Ai(S) are tabulated below 

S M M/S 2/s 

0.5 0.81 1.41 

0.4 0.82 1.51 

0.3 0.91 2.03 

0.2 1.25 3.65 

0.1 2.27 10.54 

The slow variation of M with S for S ~ 0.5 reflects the fact 
that the point of inflexion (d 2 (j>/dR 2 — 0) which defines the 
marginally Jeans-unstable fragment occurs at S = 0.5, R = 
S 2 ^ 3 Ri ~ 0.63i?i. Consequently acceleration of the collapse is 
slow at this stage. 

There is apparently an inconsistency in the above analysis, 
in the sense that - strictly speaking - a marginally Jeans-unstable 
fragment would be centrally condensed, and hence the value of J 
would be smaller than the one we have derived. However, since 
here we are mainly concerned with dynamical fragmentation, we 
do not envisage there being sufficient time for the fragment of the 
moment to relax towards detailed hydrostatic balance. Therefore it 
is appropriate to use the results we have derived on the basis of a 
uniform-density fragment. 



Appendix B. Radiative transport 
inside an optically thick fragment 

In order to estimate the temperature difference required to drive ra- 
diative energy transport from the centre of the fragment to its sur- 
face in the high optical depth limit, we introduce a dummy optical 
depth variable r' with r' — at the centre and r' = f(T) at the 
surface. For this calculation, we must allow that the temperature is 
a function of position 6(t'), and increases from 9 = T at the sur- 
face to higher values in the interior. We shall adopt the condition 
that isofhermality breaks down once a fraction / of the fragment's 
mass is heated above FT, and adopt / = 1/3, F = 2. 

Since the mean dust opacity depends on temperature (see Eq. 
Um), the effective optical-depth differential is 



(9/Tf dr 



and so the flux is given by 

4 d 
3 



For simplicity we assume that heat is deposited uniformly 
throughout the fragment, so the divergence of the flux must satisfy 



V.F 



_L A. 

r 2 dr 

l§aT f(T) 



{r 2 F(r)} 



3R 



i_ A 

^2 d7 



; { T 'V-<V)£} 



3W« 



AtyR 3 



(15) 



where we have substituted r = Rt'/t(T) and dr = Rdr' /t(T). 

Integrating this equation, and applying the boundary condi- 
tions dO/dr' = at r' = 0, and 9 = T at r' = r(T), we obtain 



3(4 - (3)H mr f{T) 



+ 



128nRaTP 



1 - 



f(T) 



(16) 



Substituting r'/f (T) = f 1/s and 9(j') < FT, gives the 
condition for approximate isothermality to break down: 



47ri? 2 crT 4 



< Tait{f,F,P) 

32 [F< 4 -« - l] 



3 [4-/?] [1-/2/3] 



(17) 



For / = 1/3 and F = 2, we obtain the values 



1 LJli 



1.00 
47.9 



1.50 
38.3 



2.00 
30.8 



Appendix C. Gas-cooling by dust 
versus gas-cooling by lines 

Consider a molecule with mass m x and moment of inertia Ix- Its ro- 
tational energy levels are at Ej = J( J+ l)h 2 /2I X ~ J 2 h 2 /2I X , 
and they will be significantly excited up to -Bj hiehes , = TkT 
(where T is a factor of order unity), so we can put 



tAishcsi 



[2ThkT] 



1/2 



hv, 



h 

The rotational lines have 

JT1 2 _ 
7-.J-1) ~ j 



Jh 

27r/x 



and hence 



hv, 



< hv 



( ^highest — * J highest ~ 1 ) 



2TkT 



To calculate the maximum possible cooling rate due to 
this molecule, we assume that all the lines up to and including 
(Jhighest — » Jhigtai - 1) are optically thick. 

Under this circumstance, the intensity in a line is given by 
the Planck Function. Again, since we are after a maximum pos- 
sible cooling rate, we can replace the Planck Function with the 
Rayleigh- Jeans approximation, B V {T) ~ 2kTv 2 /c 2 . Since the 
peak of the Planck Function is at hv ~ 5kT, whereas the highest 
frequency optically thick cooling line is at hv, , , , , — 

n J r J ° I ■'highest -> •'highest -1 ) 

2TkT j Jhighest, this is actually a fairly good approximation as long 
as Jhighea > O.4.F. 

The width of a line is determined by the velocity dispersion. 
For the thermally supported low-mass clumps with which we are 
concerned the velocity width is given by 



Av 



'2kT~ 


1/2 


" 2m " 


1/2 




= 4a 










.7rm x _ 
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If we consider a clump of mass M and radius R, then the max- 
imum cooling rate due to a single line from this particular molecule 
is 



AnR 2 7T [B„(T)-B„(T.)] — 

c 



AR 2 k [T - T ] a 



2m 



Tvm x 



1/2 



v=jh/2irl x 
W 3 



(18) 



(19) 



where T„ is the temperature of the background radiation field. 

The total cooling rate due to the molecule is obtained by sum- 
ming over all the lines: 



j=j< 



Cx 



E 

j=i 

R 2 k[T-T \a 

7T 



2m 



7rm x 



1/2 



c/x 



3 J 4 

^ highest 



4^ 2 i? 2 m 2 fc[T - 7> 5 

7TC 3 ?l/v 



2m 



7rm x 



1/2 



(20) 

(21) 
(22) 



From the results of Neufeld, Lepp & Melnick (1995), we in- 
fer that at high densities (i.e. in the optically thick limit) 12 C le O 
contributes ~ 4 % of the net cooling, so we put m x = m co ~ 
5 x 10~ 23 g, 7x = Jco 
cooling rate in the form 



6 x 10 ' gem s 1 , and write the net 



m co 



1/6 ^ 

rh 



(25) 



With J = 2.4, Z im = 0.01, Af = 25, J 7 = 2, r dust = 10~ 5 
cm, p dust = 3 g cm -3 , 7 = 5/3, and m = 4 x 10~ 24 g, this 
reduces to 



R > -Rcril 



0.3 pc 



T - T du 



T — T„ 



1/3 



(26) 



With a view to making 7?^, as small as possible, we put 
(T - T„) ~ 10K (e.g. T„ ~ 3K and T ~ 13K). Then to reduce 
Ran to 0.1 pc requires (T — T dllsl ) < 0.3K, and to reduce it to 0.03 
pc requires (T — T dust ) < 0.01K. It would be an extraordinary co- 
incidence if, in a regime where gas-cooling by lines was dominant, 
thermal balance were to deliver such a small value of (T — T dllsl ). 
We conclude that for fragments with 7? <, 0.1 pc, gas-cooling is 
dominated by dust. 
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A/"Cco 

4A/^ 2 T? 2 m 2 fc[r-r ]a 5 

lTC 3 hIco 



2m 



7rm C o 



1/2 



(23) 



where M — 25 is intended to represent the contribution from 24 
other similar molecular lines from different molecules and/or iso- 
topic variants. This net gas-cooling rate due to (molecular) lines 
should be compared with the gas-cooling rate due to dust. 

It turns out that the two cooling rates are comparable in the 
relatively low density regime where the gas-cooling rate due to dust 
is dominated by the rate of transfer of thermal energy from the gas 
to the dust, i.e. 



M r 8fcT ' 



1/2 



3pZ& 



4r dus t/5 dus t 



1 — exp 



-75K 



T 



k[T~ T dusl ] 

[7-1] 



(24) 



rather than the rate at which the dust reradiates this energy (C ad ). 
Additionally, at the low temperatures with which we are mainly 
concerned, we can simplify the analysis by setting the accommo- 
dation coefficient (the term in braces, { }) to unity. 
Putting 



feT _ 2 _ GM 

rh ~ a ~ JR ' 

the condition for line cooling to dominate over dust cooling be- 
comes 



R > 7?„ii 



9^7 2 ?i^dusi7co [T — Tdust] 
32A/"J«GV dus ,p dl , s , [ 7 - 1] [T - T B ] 



1/3 
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